Oscillatory mean-field dynamos with spherically symmetric, isotropic a 
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Until recently, the existence of oscillatory mean-field dynamos of the a^-type with spherically 
symmetric and isotropic a was an open question. We find such dynamos by means of an evolutionary 
strategy, and we illustrate the spectral properties of the corresponding dynamo operators. 
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For decades, homogeneous dynamos have been the sub- 
ject of purely theoretical research that tried to explain 
magnetic field generation in such inaccessible regions as 
the Earth's deep interior, the sun's convection zone, or 
the spiral arms of galaxies. In 1999, this situation has 
changed with the first successful dynamo experiments at 
the sodium facilities in Riga [|l], |^ and Karlsruhe ||. 
Now, dynamos run in laboratory, and their kinematic 
behavior can be predicted with an error margin of a few 
percent Dynamo theory is not unaffected by those 
developments in dynamo "engineering" . Whereas in the 
past the relevance of a dynamo model was ultimately 
judged by its applicability to real cosmic bodies, it is 
now sensible to "tailor" particular models which could 
show interesting features, e.g. field reversals or mode 
switching, in the laboratory. 

These developments give also new impetus on spec- 
tral analysis of dynamo operators. One of the simplest 
dynamo models is a mean-field dynamo with a spheri- 
cally symmetric and isotropic helical turbulence parame- 
ter a. In particular, the case with constant a is one of the 
rare models in dynamo theory that can be solved semi- 
analytically . Though the construction of a spher- 
ically symmetric, isotropic a^-dynamo seems away from and the source-free condition 
realistic cosmic dynamos with their, for instance, typical 
antisymmetry with respect to the equatorial mid-plane, 
the basic analysis of this simplest case of a^-dynamos is 
a necessary starting point for the investigation of more 
realistic dynamos. Despite its simplicity, the spherically 
symmetric, isotropic a^-dynamo exhibits a rich spectral 
structure if a is allowed to vary with the radius |8[ ^. 
In particular, complex eigenvalues have been found for 
a number of models ^ 11, [l^. However, 



radial dependence a{r) so that the eigenmode for I = 1 
has zero growth rate and non-zero frequency and that all 
other modes have negative growth rates. This task fits 
into a class of recently treated inverse problems concern- 
ing the determination of a{r) from the demand that cer- 
tain spectral properties (a few given eigenvalues, |^ ^J) 
or spatial properties of the magnetic field ("hidden dy- 
namos", [^) are fulfilled. 

In the following, we will shortly describe the forward 
solver that is used for the determination of the eigenval- 
ues, and the evolutionary strategy (ES) which we employ 
to solve the inverse problem formulated above. 

We start with the induction equation for a mean-field 
dynamo model with a spherically symmetric a-coefficient 
(for the basics of mean- field dynamos see 0). Inside a 
sphere of radius R the electrical conductivity a is con- 
stant, whereas it is zero in the exterior. In the interior 
of the sphere the magnetic field B has to satisfy the in- 
duction equation 



V = V X (aB) + AB 

at /iof 



V • B = 



(1) 



(2) 



At the boundary r = R, the magnetic field has to 
match continuously to a potential field in the exterior. 

As usual |7|, we decompose B into a poloidal and a 
toroidal part, 



it always 

turned out that other modes with higher degrees I of 
the spherical harmonics had larger growth rates than the 
considered oscillatory mode (cf. |13| ). An explicit exam- 
ple of a spherically symmetric, isotropic a^-dynamo with 
the dominating mode being oscillatory does not exist up 
to present. 

Thus motivated, our concrete goal in the present pa- 
per is to solve the following inverse problem: Find, for 
an a^-dynamo with spherically symmetric, isotropic a, a 



B = -V X (r X V5) - r X Vr , 



(3) 



with the defining scalars S and T expanded in spherical 
harmonics according to 



I — I ■m— — l 
oc I 

T{r,e,cj,) tr{r)Yr{e,c^)exp{Xit) . (5) 



t — l m— — t 



'Electronic address: 
t Electronic address: 



F.Stefani@fz-rossendorf.de 



. GerbethOfz-rossendorf . de 



For the remaining of the paper, we will measure the 
length in units of R, the time in units of iioaR^, and the 
parameter a in units of {^q(jR)~^. 
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Using Eqs. (3-5), the induction equation (1) can be 
transformed into the eigenvalue equation system: 

1 1(1 + 1) 
= ;^(^*') ^si+a{r)ti, (6) 

I d f d r \ d \ 

r ar \dr ar ) 

— {ti~a{r)si). (7) 

In our particular case of a spherically symmetric, 
isotropic a there is no coupling between field modes dif- 
fering in the degree / of the spherical harmonics. The or- 
der m of the coefficients of the defining scalars has been 
skipped as it does not show up in the equations. The 
radial wavenumber is labeled by the index n. 

The boundary conditions at r = 1 are as follows: 

^\r=l + il + l)siil)=ti{l) = 0. (8) 

ar 

A shooting technique and a fifth order Runge-Kutta 
method was used to solve this system numerically, uti- 
lizing and adapting standard routines from Numerical 
Recipes [^5). This code has been validated by compari- 
son with exact results (available only for a = constant), 
and with the results of other codes, including an integral 
equation solver []l6| . 

For the description of the inverse spectral solver we will 
be very short as it has been described in detail in j|, 1). 
We use an evolutionary strategy (ES), with a "popula- 
tion" of "individuals" (parameter vectors that describe 
the function a(r)) evolving in a "fitness" landscape ac- 
cording to principles taken from biology. Finally the evo- 
lution is stopped when the population has gathered in 
the (hopefully) global maximum of the fitness landscape. 
Note that in every step of the evolution the eigenvalue 
equation system (6-8) is solved correctly, without any 
compromise with other functionals to be minimized. 

This evolutionary strategy is very robust, and it can 
easily be used for the solution of different inverse spectral 
problems |^ . For the solution of our task we define the 
fitness F[a\ of a function a{r) according to: 

F[a] = -(5R(Ai,iH))2 

-(l + e + tanh(3(Ai,iH) -/))-! 

-(l + e + tanh(-3?(A2,i [«])))-! 

-(1 + e + tanh(-3?(A34[a])))-i . (9) 

What is the motivation for this choice of the fitness? The 
first term, to begin with, tries to keep the growth rate 
5ft(Ai^i[Q:]) of the eigenmode with {I = l\n = 1) as close as 
possible to zero. The second term makes the frequency 
9(Ai_i[a]) of this mode repel from the parameter / which 
will be considered as variable. The third and the fourth 
term are chosen so as to force the growth rates for (Z = 
2\n = 1) and {I — 3\n — 1) to be less than zero. In a 
strict sense one should also demand the growth rates for 
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FIG. 1; The functions a{r) resulting from the evolutionary 
strategy with fitness function F[a] according to Eq. (9) for 
four difi^erent values of /. 

all the higher values of I and n to be less than zero, but 
this can also be checked a-posteriori. The small value e, 
which is only used to avoid numerical overflows, is chosen 
as 10"^. 

It should be mentioned that the ES can be combined 
with a certain regularization of the functions a(r). By 
means of a regularization parameter we can adjust the 
allowed mean quadratic curvature of a{r) It turns 
out that if we start with a large regularization parameter 
that keeps the mean quadratic curvature of a{r) small 
we do not get any solution of our problem. Only with 
a rather small regularization parameter that allows more 
curvature of a{r) the ES yields solutions. 

In Fig. 1 we represent, for four different values 
of the parameter /, the functions a(r) resulting as 
solutions of the ES. Note that the functions a{r) change 
their sign at least once. As we started with large 
regularization parameters without getting solutions, the 
curves in Fig. 1 can be considered as the "simplest" 
ones (in the sense of minimal mean quadratic curva- 
ture) that fulfill our spectral demands. The analytical 
expressions for the curves are given in the following table: 
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a(r) 








0.7 


-19.88-1-347.37 




-656.71 


j' '^ ^ 


h335.52 




2.0 


-20.95-f 399.40 




-765.69 




1-387.61 


r4 


5.0 


-21.46-h426.41 




-806.73 




1-392.28 


r4 


9.0 


-23.93-h594.35 


r2 


-1243.02 




-665.27 


^4 



For these functions a{r), the growth rates of the 
eigenmodes with (/ — 1...3\n — 1) and the frequency of 
the eigenmode with {I — l\n — 1) are as follows: 



/ jR(Ai,i[a]) a(Ai,i[a]) SR(A2,iM) g?(A3,i[a]) 



0.7 


0.00 


4.98 


-5.08 


-12.02 


2.0 


0.00 


6.02 


-5.35 


-13.02 


5.0 


0.01 


7.24 


-3.64 


-10.82 


9.0 


-0.11 


9.30 


-1.42 


-9.32 
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FIG. 2: Special case / = 5.0. Growth rates for for the eigen- 
functions with {I — l\n = 1,2) and (l = 2...6|n — 1). 
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FIG. 3: Special case / = 5.0. Growth rates for the eigenfunc- 
tions with {I = l\n — 1, 2) and {I = 2...5|n = 1). Details close 
to the critical point C — 1. 



These values have been validated in an a-posteriori 
analysis of the spectrum for the functions a{r). 

Evidently, the ES delivers functions a{r) that fulfill the 
spectral demands formulated above. The {I — l\n = 1) 
eigenmode with zero growth rate is oscillatory, while the 
remaining eigenmodes have growth rates less than zero. 
To be on the safe side, we have also computed the growth 
rates for modes with / > 3, which also turned out to be 
less than zero. 

For the sake of illustration, let us concentrate in the 
following on the spectral properties for the special case 
with / = 5. In Fig. ^ we show the growth rates for the 
eigenmodes with (l — l\n = 1,2) and {I = 2...6\n — 1) 
as functions of C, which is used to scale the magnitude 
of a{r) (for C = we have the free decay case, C = 1 
corresponds to the obtained functions a{r)). Details of 
this plot close to the critical point C = 1 are shown in 
Fig I 

At C = 0, all modes start as non-oscillatory modes 
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FIG. 4: Special case / = 5.0. Growth rates for the eigenfunc- 
tions with (l = l|n — 1...7). 
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FIG. 5: Special case / — 5.0. Frequencies for the eigenfunc- 
tions with {I = l|n = 1...7). 



(with purely real eigenvalues). At C = 0.818, the two 
modes with (l = l\n = 1) and {I ^ l\n ~ 2) merge and 
continue as a pair with complex conjugate eigenvalues 
(only a single line is shown). Interestingly enough, at 
C = 1.097 this pair splits off again and the two modes 
with {I = l\n = 1) and {I = l\n = 2) continue separately, 
again with purely real eigenvalues. In between, at C — 
1.00 the pair crosses the abscissa, hence dynamo action 
occurs. All modes with higher I and n have still negative 
growth rates at C = 1.00. 

In Figs. ^ and || the growth rates and frequencies for 
the modes with {I — l\n — 1...7) are plotted. These fig- 
ures give a flavor of the complexity of the spectrum, with 
its merging and splitting points of modes with different 
n. 

In this respect, another aspect should be addressed. 
Whereas the spectra for the functions a{r) for / = 
0.7, 2.0, 5.0 are rather similar, there is a transition if we 
go to / = 9.0. Figure Hshows the corresponding spec- 
trum. Whereas in Fig.T3 we observed a merging of the 
modes with n = 1 and n = 2, we get now a merging 
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FIG. 6: Special case / = 9.0. Growth rates for the eigenfunc- 
tions with {I = l|n — 1...7). 
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FIG. 7: Modified functions a{r) of the type obtained for 
/ = 5, but with varying coefficient of the term, that are 
oscillatory dynamos. 



of the modes with n = 2 and n — 3. For larger C, the 



modes with n = 1 and n = 4 come together, meeting 
at the same point with the common line of the complex 
conjugated n = 2(3) modes. 

Evidently, the dynamo operators connected with the 
obtained functions a(r) show very interesting spectral 
features which warrants further analysis. 

A complete characterization of the a(r)-profiles that 
fulfill our spectral demands would require an extensive 
parameter study. At least we can get a certain feeling 
on the broadness of the "corridor" for a limited type of 
deformations. For this purpose, we start with the func- 
tion a{r) ^ -21.46 + 426.41r2 - 806.73r3 + 392.28r4 
which was the solution for / = 5. Now let us change the 
coefficient of the term, leaving all other coefficients 
unchanged. Figure ^ shows the upper and the lower lim- 
its of the corresponding deformations that still belong to 
our oscillatory class of dynamos. Above the upper limit, 
the {I — 2\n — 1) mode starts to dominates. Below the 
lower limit, the {I — 4|ri — 1) mode dominates. 

In summary, we have found a class of oscillatory mean- 
field dynamos with spherically symmetric, isotropic a 
working in a finite volume with homogeneous conductiv- 
ity. The obtained functions a(r) are smooth and by no 
means exotic or of a pronounced layer type |ll| , but 
they are characterized by at least one change of sign along 
the radius. In the sense of minimal averaged quadratic 
curvature, our solutions seem to be the simplest ones. 

As for the back-reaction regime, an interesting be- 
haviour of the considered dynamo models is conceivable. 
Even a slight modification of the functions a(r) due to 
the action of the Lorentz forces could trigger a transi- 
tion to a non-oscillatory mode (either with ^ = 1 or with 
higher I). 

We do not claim any particular astrophysical relevance 
of our result. However, keeping in mind that the Karls- 
ruhe dynamo experiment is an a^-dynamo (although 
with an anisotropic a-tensor) one could imagine a gener- 
alization of our method to similar laboratory dynamos. 
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